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1 . 1  BACKGROUND 

The  U.S.  Coast  Guard  conducted  a  comprehensive  data 
collection  effort  from  May  1979  to  May  1980  to  quantify  an 
apparent  temporal  instability  in  the  St.  Marys  River  Loran-C 
Time  Difference  (TD)  grid.  Detailed  analyses  of  the  resulting 
Loran-C  data  base  are  documented  in  this  report. 

TASC ' s  involvement  with  the  St.  Marys  River  Loran-C 
mini-chain  began  with  the  calibration  of  a  spatial  TD  grid 
prediction  model  based  on  data  collected  during  September  and 
October  1977  (Refs.  1  and  2).  A  limited  amount  of  data  collected 
monthly  between  November  1977  and  March  1978  exhibited  temporal 
TD  variations  as  large  as  0.4  psec  peak-to-peak  at  some  data 
collection  sites.  Because  the  utility  of  the  TD  grid  prediction 
model  and  the  mini-chain  depends  on  the  assumption  of  a  stable 
grid,  it  was  important  that  the  validity  of  the  observed  temporal 
variations  be  confirmed.  In  particular,  collection  of  additional 
data  under  more  carefully  controlled  conditions  (e.g.,  well- 
sheltered  receivers  and  ideal  grounding)  was  judged  to  be  necessary. 
The  need  for  confirmation  was  further  highlighted  by  a  theoretical 
sensitivity  analysis,  which  revealed  that  the  observed  TD  varia¬ 
tions  could  not  be  explained  by  expected  seasonal  variations 
in  signal  propagation  conditions  (Ref.  3).  In  order  to  resolve 
the  discrepancy  between  the  theory  and  observations,  it  was 
recommended  that  the  U.S.  Coast  Guard  initiate  the  yearlong 
data  collection  effort  outlined  in  Ref.  4.  Loran-C  Data-Base 
Management  Software  was  subsequently  designed  and  implemented 
at  TASC,  to  support  the  editing  and  analysis  of  the  collected 
data  (Ref.  5 ) . 
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The  data  consist  of  TD  samples  and  associated  signal- 
quality  indicators,  recorded  nominally  every  15  min  at  three 
fixed  sites:  DeTour  Village,  Dunbar  Experimental  Forest,  and 

•Ar 

Pt.  Iroquois  (see  Fig.  1.1-1).  Two  Internav  LC-204  Loran-C 
receivers  are  located  at  both  the  DeTour  and  Iroquois  sites, 
and  one  Internav  LC-204  and  one  Magnavox  AN/BRN-5  receiver  are 
located  at  the  Dunbar  site.^  The  BRN-5  receiver  is  interfaced 
with  a  cesium  frequency  standard  to  also  provide  measurements 
of  Time-of-Arrival  (TOA) .  The  following  additional  data  are 
collected  to  support  the  analyses:  System  Area  Monitor  (SAM) 
TDs  from  the  Austron-5000  receiver  used  for  chain  control, 
Local  Phase  Adjustments  (LPAs)  from  the  Calculator  Assisted 
Loran  Controller  (CALOC)  system,  and  meteorological  data  from 
the  National  Weather  Service  Station  at  Sault  Sainte  Marie, 
Michigan. 


In  addition  to  the  15-min  TD  data,  which  are  termed 
high-density  data,  the  U.S.  Coast  Guard  collects  low-density 
data  for  in-house  analyses.  The  low-density  data  are  recorded 
with  a  microprocessor-based  system  which  is  independent  of  the 
Texas  Instruments  Silent  700  cassette  tape  recorder  used  to 
obtain  the  high-density  data.  The  low-density  data  consist  of 
five  15-min  samples  recorded  twice  daily,  which  are  then  averaged 
and  transmitted  over  a  telephone  data  link  in  a  weekly  message. 

A  similar  low-density  data-recording  system  is  under  development 
by  the  U.S.  Coast  Guard  for  use  in  harbor  and  harbor  entrance 
Loran-C  navigation  experiments. 


*Referred  to  as  DeTour,  Dunbar,  and  Iroquois  in  this  report. 
tReferred  to  as  LC-204  and  BRN-5  in  this  report. 
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Objectives  and  Technical  Approach  for 
St.  Marys  River  Temporal  Loran-C 
Data  Analyses 


Evaluate  the  low-density  data  analysis 
approach  to  establish  an  efficient  Loran-C 
monitoring  strategy  for  harbors. 


To  meet  these  objectives,  the  raw  Loran-C  data  are  first 
organized  and  edited,  using  the  TASC  Loran-C  Data-Base  Manage¬ 
ment  Software.  Editing  marks  are  set  by  the  software  to 
indicate  the  quality  of  each  TD  sample,  thus  making  it  straight 


L 


forward  to  extract  "good"  and/or  "bad"  data  without  affecting 
the  master  data  base.  Specifically,  a  criterion  for  editing 
is  selected  for  the  high-density  analyses,  and  all  correspond¬ 
ing  data  are  extracted  and  analyzed  (see  Fig.  1.2-1).  Availa¬ 
bility  of  the  high-density  analysis  results  then  permits  the 
low-density  evaluation  to  be  limited  to  a  subset  of  all  possible 
site/receiver/TD  combinations.  Both  edited  and  unedited  TD 
data  are  extracted  for  the  selected  combinations  to  determine 
the  effect  of  outliers  on  the  low-density  approach. 

It  should  be  kept  in  mind  that  the  primary  objective 
of  the  high-density  data  analyses  is  to  quantify  TD  grid  insta¬ 
bility.  Although  the  determination  and  modeling  of  cause/ 
effect  relationships  is  beyond  the  scope  of  the  current  effort, 
the  analyses  are  selected  and  presented  in  a  manner  which  facili¬ 
tates  possible  future  studies  by  the  U.S.  Coast  Guard.  A  compari¬ 
son  is  made  between  the  data  analysis  results  herein  and  the 
theoretical  sensitivity  analysis  results  from  Ref.  3,  to  provide 
a  basis  for  future  research. 


1.3  REPORT  OVERVIEW 

A  description  of  the  Loran-C  data  base  and  editing 
procedure  is  provided  in  Chapter  2.  The  high-density  data 
analysis  results  are  given  in  Chapter  3  and  the  low-density 
data  analysis  approach  is  evaluated  in  Chapter  4.  The  study 
conclusions  are  contained  in  Chapter  5,  together  with  a 
recommended  sarapling/editing/analysis  strategy  for  harbor 
Loran-C  monitor  data. 

The  TASC  Loran-C  data  base  is  being  delivered  to  the 
U.S.  Coast  Guard  on  three  magnetic  tapes  which  supplement  this 
report.  The  tape  layouts  are  described  in  Section  2.6  and  the 
record  formats  are  given  in  the  appendices. 
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DESCRIPTION  OF  DATA  BASE  AND  EDITING  PROCEDURE 
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2 . 1  INTRODUCTION 

Data  from  the  St.  Marys  River  temporal  TD  instability 
experiment  are  provided  to  TASC  on  11  half-inch  magnetic  tapes 
with  the  record  format  listed  in  Appendix  A.  The  tapes  contain 
a  total  of  137  site  files,  286  SAM  files,  and  11  weather  files, 
which  are  sorted  and  ordered  by  the  Loran-C  Data-Base  Management 
Software.  The  record  formats  are  changed  to  those  presented 
in  Appendices  B  and  C,  which  include  six  editing  mark  characters 
for  every  TD  and  TOA  sample. 


2 . 2  EDITING  PROCEDURE 

The  editing  marks  indicate  the  quality  of  the  data, 
each  mark  corresponding  to  a  different  editing  operation  as 
indicated  in  Table  2.2-1.  The  manner  in  which  the  editing 
marks  are  determined  is  detailed  in  Ref.  5  and  summarized  below. 


2.2.1  Editing  Mark  1:  Loran-C  Transmitter 


Editing  Mark  1  indicates  which  Loran-C  transmitter  is 
associated  with  the  TOA,  or  which  secondary  transmitter  is 
associated  with  the  TD.  In  the  cases  of  the  LC-204  and  BRN-5 
receivers,  the  transmitter  is  determined  by  examining  the  first 
two  digits  of  the  TD.  If  the  digits  do  not  match  one  of  the 
three  possibilities  --  "11"  for  TDX,  "22"  for  TDY,  or  "33"  for 
TDZ  --  Editing  Marks  1  to  6  are  set  to  "?". 
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TABLE  2.2-1 

INTERPRETATION  OF  EDITING  MARK  VALUES 


EDITING 

MARK 

OPERATION 

VALUES* 

INTERPRETATION 

1 

Loran-C  Transmitter 

M 

TO  AM 

X 

TOAX  or  TDX 

Y 

TOAY  or  TDY 

z 

TOAZ  or  TDZ 

Signal  Quality 

0 

Good  Data  Quality 

Check 

1 

Control  Unit  Switch  on  B 

2 

Insufficient  Signal 

3 

Both  Characteristics  1  and  2  Apply 

1 

Tolerance  Check 

0 

Data  Differs  froa  Seasonal  Mean 

by  <  0.5  Msec 

1 

0.5  usee  <  Difference  <_  1.0  p»ec 

2 

Difference  >  1.0  usee 

u 

Manual  Editing 

0 

No  Manual  Editing 

6 

Duplicate  Data 

5 

Differential 

Differential  Correction  Computed 

Correction 

0 

No  Interpolation  Required 

Computation 

A 

Interpolation  Interval  <_  6  hr 

B 

6  hr  <  Interval  <_  24  hr 

C 

Interval  >  24  hr 

Differential  Correction  Not  Computed 

l 

Uncertain  LPA  Times 

2 

Bad  SAM  TD  Data  Based  on  Editing 

Marks  1  to  4  and  6 

3 

No  LPA  and/or  SAM  TD  Data  Available 

6 

Outlier  Editing 

0 

Data  Is  Mot  an  Outlier 

1 .2 ,3 .4 

Data  Differs  from  Daily  Trend  Line  by 

>0.05,  0.10,  0.15,  0.20  nsec  (Second 

Iteration) 

5,6 

Data  Differs  from  Trend  Line  by  >0.15, 

0.20  psec  (First  Iteration; 

S 

Less  Than  20  Samples  per  Day  Available  to 

Compute  Trend  Line 

U 

Bad  Data  Based  on  Editing  Marks  1  to  4 
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2.2.2  Editing  Mark  2:  Signal-Quality  Check 


Editing  Mark  2  summarizes  various  Loran-C  signal- 
quality  indicators  which  appear  in  each  data  record.  The  in¬ 
dicators  which  are  checked  are  SWITCH  and  SUFFSIG  for  the  LC-204, 
SUFS1G  for  the  BRN-5,  and  MODE  for  the  Austron-5000  receiver 
(see  Appendix  A).  Editing  Mark  2  is  the  only  editing  mark 
which  is  not  based  on  the  TD/TOA  data  itself. 

2.2.3  Editing  Mark  3:  Tolerance  Check 

Editing  Mark  3  indicates  whether  or  not  the  TD  lies 
between  expected  lower  and  upper  bounds,  which  are  selected 
based  on  U.S.  Coast  Guard  low-density  data  available  at  the 
beginning  of  the  study.  Specifically ,  the  TD  is  compared  to 
the  mean  of  the  low-density  dafa  for  the  period  May  1979  to 

*iAr 

February  1980  (Ref.  5).  Three  tolerance  bands  (see  Table 
2.2-1)  are  distinguished  to  permit  the  editing  of  gross  out¬ 
liers,  while  retaining  the  capability  to  accommodate  TD  varia¬ 
tions  during  the  March  1980  to  May  1980  time  period.  Editing 
Mark  3  serves  primarily  to  identify  TDs  which  are  affected  by 
receiver  cycle  jumps  and  transmitter  and  control  malfunctions. 

2.2.4  Editing  Mark  4:  Manual  Editing 

Editing  Mark  4  is  employed  to  manually  flag  suspect 
data  which  are  not  flagged  automatically  by  other  editing  marks. 
The  large  volume  of  data  makes  it  impractical  to  employ  the 
Test  Director's  Log  for  manual  editing.  Consequently,  the 
only  use  made  of  Editing  Mark  4  is  to  identify  duplicate  data 
which  arise  either  in  the  transferral  of  tape  cassette  data  to 
half-inch  magnetic  tape  by  the  U.S.  Coast  Guard,  or  in  the 


*The  mean  of  the  high-density  TD  data  for  the  entire  year 
differs  from  this  value  by  less  than  0.05  psec. 
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processing  of  the  half-inch  magnetic  tapes  by  TASC.  The  absence 
of  further  manual  editing  does  not  detract  from  the  data  quality 
in  view  of  the  comprehensive  automatic  editing  employed. 

2.2.5  Editing  Mark  5:  Differential-Correction  Indicator 

Differential  corrections  are  computed  and  included  in 
the  site  data  records  to  show  the  difference  between  measured 
SAM  TDs  and  the  Controlling  Standard  TD  (the  reference  used 
for  chain  control).  Differential  corrections  can  be  added  to 
the  measured  site  TDs  in  post-time  to  simulate  perfect  chain 
control.  When  site  and  SAM  TDs  are  not  measured  at  the  same 
times,  the  SAM  TDs  are  estimated  by  linear  interpolation  (Ref. 
5).  Care  is  taken  to  include  the  effects  of  LPAs  and  receiver 
averaging  times  (~50  sec)  in  the  interpolated  TDs.  An  offset 
of  850  sec  in  the  LC-204  recording  times  relative  to  the  actual 
sampling  times  is  also  accounted  for  in  this  procedure. 

Editing  Mark  5  indicates  the  length  of  the  interpola¬ 
tion  interval  used  in  the  differential-correction  computation 
(see  Table  2.2-1).  Intervals  greater  than  a  quarter  of  the 
diurnal  cycle  (6  hr)  are  judged  to  be  excessive,  because  the 
linear  approximation  no  longer  holds.  Table  2.2-1  indicates 
circumstances  under  which  no  differential  correction  can  be 
computed.  Noteworthy  in  this  regard  is  the  period  from  22 
January  1980  to  31  March  1980,  for  which  LPAs  were  provided  to 
TASC  in  hand-written  form  with  time  to  the  nearest  minute, 
rather  than  on  the  CALOC  system  printout  with  time  to  the 
nearest  second.  Editing  Mark  5  is  set  to  "3"  in  this  case 
because  the  accuracy  of  the  LPA  data  has  not  been  verified  by 
TASC . 


2-4 


V, 


2.2.6  Editing  Mark  6:  Outlier  Editing 

Editing  Mark  6  is  employed  to  flag  data  which  are 
inconsistent  with  the  short-term  behavior  of  the  TD  time  series, 
but  are  not  flagged  by  other  editing  marks.  These  outlying 
data  represent  isolated  "events"  caused  by  low  signal-to-noise 
ratio,  transraitter/receiver/control  malfunctions,  or  anomalous 
propagation  conditions.  Although  the  outliers  must  be  edited 
to  permit  an  orderly  analysis  of  typical  TD  variations,  they 
are  marked  in  the  data  base  for  use  in  possible  future  cause/ 
effect  studies. 

Outliers  are  detected  automatically  by  comparing  each 
TD  sample  to  a  trend  line  fit  to  all  data  in  a  daily  time  window 
(Ref.  5).  Samples  which  differ  from  the  trend  line  by  greater 
than  0.15  psec  are  marked  "5"  or  "6"  as  indicated  in  Table 
2.2-1,  and  a  new  trend  line  is  computed  from  the  remaining  samples. 
The  samples  are  then  compared  to  the  refined  trend  line  and 
marked  "0"  to  "4"  to  indicate  which  0.05-psec  error-bin  they 
fall  into  (see  Table  2.2-1).  This  procedure  produces  a  histo¬ 
gram  of  TD  errors  relative  to  the  trend  line. 


2.3  TD  DATA  INVENTORY 


The  operations  involved  in  editing  the  TD  data  base 
for  high-density  analyses  are  summarized  in  Table  2.3-1.  The 
editing  mark  values  which  define  the  data  base  are  noted  below: 


•  Editing  Mark  1  =  "X",  "Y" ,  or  "Z":  two 
TD  time  series  are  extracted  for  each 

site  receiver  and  three  for  the  SAM  receiver 

•  Editing  Mark  2  =  "0":  87.7  to  98.4  percent 

of  the  data  are  associated  with  good 
signal  quality 
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•  Editing  Mark  3  =  "0":  93.0  to  99.2  per¬ 

cent  of  the  data  are  within  1.0  psec  of 
the  seasonal  mean,  except  for  TDY  for 
DeTour/  LC-204/1*  which  exhibits  numerous 
cycle  jumps 

•  Editing  Mark  4  =  "0" :  duplicate  data 
are  ignored  in  the  "Total  Number  of 
Samples"  column  in  Table  2.3-1 

•  Editing  Mark  5  =  any  value:  differential 
corrections  are  not  applied  because  of 
their  small  magnitude  relative  to  site 

TD  variations  (see  Section  3.2)  and  because 
they  are  not  available  for  the  entire 
year 

•  Editing  Mark  6  =  "0":  89.2  to  98.1  percent 

of  the  data  are  within  0.05  psec  of  the 
daily  trend  lines  used  for  outlier  detection. 


The  number  of  samples  retained  after  all  editing  operations 
ranges  from  10,000  to  35,000  for  the  15  TD  time  series  (see 
Table  2.3-1). 


i 


The  quality  of  the  edited  data  base  is  judged  to  be 
excellent .  Not  only  do  all  edited  data  lie  within  a  ±0.05  psec 
band  centered  on  the  daily  trend  line,  but  the  standard  devia¬ 
tions  of  the  TD/trend  line  differences  are  between  0.015  psec 
and  0.031  psec  (see  Table  2.3-1).  The  short-term  fluctuations 
are  caused  by  atmospheric  noise  and  the  diurnal  TD  cycle,  and 
are  consistent  with  those  measured  in  the  laboratory  (Ref.  7). 

An  inventory  of  the  edited  TD  data  is  presented  in 
Figs.  2.3-1  to  2.3-7,  where  the  bars  indicate  the  number  of 
samples  per  day  (nominally,  96).  The  large  gaps  which  appear 


*Each  receiver  is  uniquely  identified  by  site/receiver  type/ 
receiver  number. 
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Figure  2.3-4  TD  Data  Inventory  for  Dunbar/LC-204/2 
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Figure  2.3-6  TD  Data  Inventory  for  Iroquois/LC-204/2 
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in  many  time  series  are  associated  with  missing  cassette  tapes, 
rather  than  edited  data.  The  U.S.  Coast  Guard  low-density 
data  can  be  used  to  fill  these  gaps  and  provide  additional 
low-frequency  (seasonal)  detail.  However,  inclusion  of  low- 
density  data  is  beyond  the  scope  of  the  current  effort  and 
would  not  significantly  affect  the  analysis  results. 


2.4  TO A  DATA  BASE 

Master,  Secondary  X,  and  Secondary  Y  TOA  data  are 
measured  with  the  BRN-5  receiver  at  the  Dunbar  site.  An  in¬ 
spection  of  the  TOA  data  shows  that  the  clock  was  reset  peri¬ 
odically  to  keep  the  master  TOA  near  a  value  of  1000  psec. 

The  numerous  TOA  jumps  introduced  by  this  procedure  make  it 
infeasible  to  remove  clock  effects  from  the  data.  Identifica¬ 
tion  of  the  jumps  is  further  complicated  by  large  gaps  in  the 
TOA  time  series,  coincident  with  the  TD  time  series  gaps  shown 
in  Fig.  2.3-3. 

Because  clock  effects  cannot  be  removed,  the  TOA  data 
are  not  analyzed  in  the  current  study.  TOA  Editing  Marks  1, 

2,  and  4  are  set  by  the  Data-Base  Management  Software,  but  the 
remaining  editing  marks  require  prior  compensation  for  clock 
phase  variations. 


2.5  LPA  DATA  INVENTORY 

LPA  data  are  manually  entered  into  the  data  base  by 
TASC ,  from  the  CALOC  system  printout.  The  number  of  LPA  samples 
per  day  for  each  secondary  transmitter  is  shown  in  Fig.  2.5-1. 
Approximately  10  LPAs  per  day,  on  the  average,  are  required  to 
satisfy  the  0.015-psec  control  tolerance. 
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2.6  DATA* BASE  TAPES  DELIVERED  TO  U.S.  COAST  GUARD 

The  46-megabyte  St.  Marys  River  Loran-C  data  base, 
with  editing  marks  included,  is  delivered  to  the  U.S.  Coast 
Guard  on  three  half-inch  magnetic  tapes  (9-track,  800  BPI , 

ASCII  code).  The  tapes  are  arranged  as  shown  in  Table  2.6-1, 
with  one  file  for  each  receiver  except  the  BRN-5,  which  requires 
a  separate  file  for  secondary  and  master  records.  The  block 
size  in  each  receiver  file  is  8190,  which  includes  63  records 
with  130  characters  each.  The  block  size  for  the  LPA  file  is 
8160,  which  includes  204  records  with  40  characters  each.  The 
record  formats  are  given  in  Appendices  B  and  C. 


TAPE 

NUMBER 


TABLE  2.6-1 

SUMMARY  OF  DATA-BASE  TAPES  u 
DELIVERED  TO  U.S.  COAST  GUARD' 


FILE 

NUMBER 


SITE 

NUMBER! 


RECEIVER 

TYPE 


BRN5  SEC 

BRN5  MAS 

LC204 

LC204 

BRN5  SEC 

BRN5  MAS 

LC204 

LC204 

LC204 

LPAs 

AUSTRON 


RECEIVER 

NUMBER 


NUMBER 
OF  DATA 
RECORDS 


2,589 
1 ,292 
19,579 
20,629 
48,017 
23,251 
20,495 
28,936 
30,040 
9,517 
151,179 


*This  table  is  repeated  on  Tape  1/File  1. 

tl  =  DeTour,  2  =  Dunbar,  3  =  Iroquois,  and  8  =  SAM. 
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QUANTIFICATION  OF  ST.  MARYS  RIVER 
LORAN-C  TD  GRID  INSTABILITY 


3.1  INTRODUCTION 

The  edited  data  base  described  in  Section  2.3  is 
employed  in  this  chapter  to  quantify  TD  grid  instability  in 
the  St.  Marys  River  Loran-C  mini-chain.  Quantification  is 
based  on  the  following  non-parametric  descriptions  of  the  data: 


TD  time  series  plots  and  related  statis¬ 
tics  to  exhibit  seasonal,  monthly,  and 
diurnal  variations 

Autocorrelation  and  cross-correlation 
functions  to  identify  TD/TD,  TD/tempera- 
ture ,  and  TD/refractive  index  relationships 

Power  spectral  densities  to  examine  the 
frequency  content  of  the  TD  variations 

LPA  statistics  to  quantify  SAM  TD  varia¬ 
tions  in  the  case  of  no  chain  control 

Comparison  of  theoretical  and  measured 
TD  variations  to  provide  a  basis  for 
possible  future  studies. 


The  high-density  analyses  parallel  the  low-density  analyses 
conducted  by  the  U.S.  Coast  Guard,  where  possible. 


3.2  TD  TIME  SERIES  PLOTS  AND  RELATED  STATISTICS 


3.2.1  Seasonal  Variations 


Seasonal  TD  variations  are  quantified  by  the  mean  and 
standard  deviation  of  the  data  in  3-day  time  windows  (see 


Figs.  3.2-1  to  3.2-4).  In  this  study,  seasonal  variations  are 
defined  as  the  variation  in  the  3-day  mean  over  the  data  collec¬ 
tion  period.  The  standard  deviation  gives  an  indication  of 
shorter-term  TD  variations  relative  to  the  seasonal  variations. 

To  avoid  misinterpretation  of  the  results,  no  point  is  plotted 
for  3-day  windows  which  contain  fewer  than  20  samples.  The  TD 
is  shown  relative  to  the  average  value  of  the  time  series, 
i.e.,  the  mean  of  the  3-day  means  (see  Table  3.2-1).  The  three 
redundant  TDs  --  TDY  for  DeTour/LC-204/1 ,  TDY  for  Dunbar/BRN-5/1 , 
and  TDZ  for  Iroquois/LC-204/1  --  are  not  included  here,  but 
are  discussed  in  Section  3.3.2. 


An  examination  of  Figs.  3.2-1  to  3.2-4  leads  to  the 
following  conclusions: 


•  Seasonal  TD  variations  are  less  than  0.4 
psec  peak-to-peak  (p-p)  at  the  DeTour, 
Dunbar,  and  Iroquois  sites,  consistent 
with  previous  observations  (Refs.  2  and  3) 

•  TD  values  in  May  1979  and  May  1980  are 
approximately  equal,  suggesting  that  the 
seasonal  meteorological  cycle  causes  the 
variations 

•  TD  variations  during  monthly  periods  are 
larger  in  the  winter  than  in  the  summer, 
consistent  with  signal  propagation  theory 
(Refs.  3  and  9) 

•  The  3-day  standard  deviations  (shown  by 
dotted  curves  in  Figs.  3.2-1  to  3.2-4) 
are  less  than  0.025  psec,  which  is  small 
relative  to  the  variation  of  the  3-day 
means 

•  The  standard  deviation  of  the  3-day  means 
(see  Table  3.2-1)  is  larger  for  TDY  than 
for  TDX  and  TDZ,  i.e.,  0.069  psec  average 
compared  to  0.032  psec  average 

•  SAM  TD  variations  are  less  than  0.02  psec 
p-p  except  for  TDY  in  February  (see  Fig. 
3.2-4),  indicating  that  chain  control 
tolerances  are  satisfied. 
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Figure  3.2-1  TDX  Seasonal  Time  Series 
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Figure  3.2-2  TDY  Seasonal  Time  Series 
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FFERENCE  (^sec)  TIME  DIFFERENCE  (^sec)  TIME  DIFFERENCE 


TIME  DIFFERENCE  (^sec)  TIME  DIFFERENCE  (nsec)  TIME  DIFFERENCE  (^sec) 


3-DAY  MEAN  t  1  a 
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TIME  (month) 
a)  Time  Difference  X 


TABLE  3.2-1 

SEASONAL  TD  STATISTICS 


i 


RECEIVER 

STATISTICS 

OF  3 -DAY  TD 

MEANS  (Msec) 

SITE 

TYPE/ 

NUMBER 

1 1  nt 

DIFFERENCE 

MEAN 

STANDARD 

DEVIATION 

PEAK-TO-PEAK 

X 

11328.261 

0.032 

0.174 

DeTour 

LC-204/1 

Y 

22210.351 

0.106 

0.532 

Y 

22210.346 

0.082 

0.300 

DeTour 

LC-204/2 

Z 

33416.693 

0.013 

0.061 

X 

11293.276 

0.027 

0.119 

Dunbar 

BRN-5/1 

Y 

22363.361 

0.045 

0.  180 

Y 

22363.369 

0.032 

0.  189 

Dunbar 

LC-204/2 

Z 

33290.870 

0.033 

0.204 

X 

11325.374 

0.037 

0.146 

Iroquois 

LC-204/1 

z 

33156.894 

0.044 

0.235 

Y 

22413.893 

0.052 

0.232 

I roquois 

LC-204/2 

Z 

33156.915 

0.034 

0.162 

X 

11259.695 

0.002 

0.013 

SAM 

Austron/ 1 

Y 

22363.866 

0 . 005 

0.037 

Z 

33200.397 

0.003 

--  —1 

0.015 

The  seasonal  component  of  TD  variations  dominates  the  monthly 
and  diurnal  components,  as  expected  for  Loran-C  groundwave 
propagation . 

3.2.2  Monthly  Variations 

The  distinct  difference  between  summer  and  winter  TD 
variations  is  analyzed  by  focusing  on  the  months  of  July  and 
January.  Monthly  TD  variations  are  quantified  by  the  mean  and 
standard  deviation  of  the  data  in  6-hr  time  windows  (see  Figs. 
3.2-5  to  3.2-8  for  July  and  Figs.  3.2-9  to  3.2-12  for  January). 
No  point  is  plotted  for  6-hr  windows  which  contain  fewer  than 
10  samples. 

An  examination  of  the  monthly  time  series  leads  to 
the  following  conclusions: 

•  TD  variations  are  less  than  0.10  psec  p-p 
for  July,  but  as  large  as  0.20  psec  p-p 
for  January 

•  TDs  often  change  very  abruptly  in  January, 

--  e.g.,  0.10  psec  during  a  one-  or  two-day 
period 

•  Diurnal  variations  are  smaller  in  July 
than  in  January,  but  represent  a  greater 
proportion  of  the  overall  monthly  variation 
in  July 

•  The  6-hr  standard  deviations  are  small 
relative  to  the  variation  of  the  6-hr 
means 

•  The  chain  is  controlled  more  effectively 

by  the  SAM  in  January  than  in  July  (compare 
Figs.  3.2-8  and  3.2-12). 


Confirmation  of  the  expected  differences  between  summer  and 
winter  TD  variations  is  an  important  contribution  of  the  St. 
Marys  River  Loran-C  data  collection/analysis  effort. 
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Figure  3.2-7  TDZ  Monthly  Time  Series  for  July 
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TDX  Monthly  Time  Series  for  January 
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Figure  3.2-10  TDY  Monthly  Time  Series  for  January 
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Figure  3.2-11  TDZ  Monthly  Time  Series  for  January 
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3.2.3  Diurnal  Variations 


Although  diurnal  variations  are  small  compared  to 
longer-term  TD  variations,  their  magnitude  is  important  in 
selecting  pattern  monitor  sampling  rates  and  user  update  rates 
for  differential  Loran-C.  The  diurnal  TD  variations  for  July 
and  January  are  analyzed  separately  herein,  using  the  following 
procedure : 


•  Sort  the  data  by  day  and  subtract  the 
daily  mean  to  form  TD  residuals,  thus 
removing  the  day-to-day  variation 

% 

•  Sort  the  TD  residuals  by  time-of-day,  $ 

irrespective  of  day-of-month ,  and  compute  i 

the  mean  and  standard  deviation  of  the  3 

residuals  in  2-hr  time  windows  (the  mean 

is  the  diurnal  cycle)  - 

•  Construct  a  histogram  of  the  differences 
between  the  residuals  and  the  diurnal 
cycle. 


The  resulting  diurnal  cycles  and  histograms  are  presented  in 
Figs.  3.2-13  to  3.2-15  for  July  and  in  Figs.  3.2-16  to  3.2-18 

JU 

for  January. 

An  examination  of  these  results  leads  to  the  following 
conclusions : 

•  The  diurnal  cycles  are  small,  generally 
less  than  0.02  psec  p-p 


*The  time  shown  is  Greenwich  Mean  Time  (GMT),  where  0  hours 
GMT  is  7  p.m.  Sault  Sainte  Marie  local  time. 


i 
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TDX  Diurnal  Cycles  and  Histograms  for  January 
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Figure  3.2-17  TDY  Diurnal  Cycles  and  Histograms  for  January 
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•  The  diurnal  cycle  is  obscured  by  random 
variations  (i.e.,  a  large  standard  devia¬ 
tion)  for  many  time  series,  but  is  very 
distinct  for  others  (e.g.,  see  Fig.  3.2-l4a) 

•  Some  diurnal  cycles  exhibit  a  bimodal 
behavior,  indicating  a  12-hr  period 
(e.g.,  see  Fig.  3.2-17a) 

•  The  2-hr  standard  deviations  are  larger 
for  January  than  for  July,  which  is  incon¬ 
sistent  with  the  lower  atmospheric  noise 
level  known  to  occur  in  January  (Ref.  10) 

•  The  2-hr  standard  deviations  for  July 
are  generally  larger  at  night  (0  to  12 
hours  GMT)  than  during  the  day,  which  is 
inconsistent  with  the  lower  noise  levels 
known  to  occur  at  night  (Ref.  10) 

•  The  random  variations  around  the  diurnal 
cycle  can  be  characterized  by  a  normal 
distribution  with  a  standard  deviation 
of  approximately  0.02  psec. 


The  high-frequency  detail  afforded  by  the  15-min  sampling  inter¬ 
val  makes  the  St.  Marys  River  Loran-C  data  base  an  important 
resource  for  future  cause/effect  analyses  of  diurnal  variations. 


3.3  CORRELATION  AND  SPECTRAL  ANALYSES 

3.3.1  Autocorrelation  Functions  and  Power 
Spectral  Densities 

The  relative  contribution  of  seasonal,  monthly,  and 
diurnal  TD  variations  is  given  by  the  Power  Spectral  Density 
(PSD),  which  is  discussed  in  detail  in  Refs.  11  to  13.  The 
following  procedure  is  used  to  compute  the  PSD  in  the  present 
study : 
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•  Construct  a  smoothed  time  series,  consist¬ 
ing  of  2-hr  averages  of  the  original 
data,  to  reduce  computer-processing 
requirements 

•  Synthesize  TD  data  by  linear  interpola¬ 
tion  to  fill  any  time  series  gaps 

•  Remove  the  mean  and  trend  (i.e.,  a  straight- 
line  fit)  from  the  time  series  to  minimize 

JL. 

non-s tationarity 

•  Compute  the  Autocorrelation  Function 
(ACF)  of  the  time  series  by  the  technique 
described  in  Ref.  13,  which  involves  the 
Fast  Fourier  Transform 

•  Multiply  the  ACF  by  a  Kaiser  window 

(Ref.  14),  which  has  the  effect  of  smooth¬ 
ing  the  PSD 

•  Compute  the  PSD  as  the  Fast  Fourier  Trans¬ 
form  of  the  ACF. 


This  procedure  is  adopted  based  on  previous  TASC  experience  in 
data  processing,  after  conducting  various  tests  on  the  present 
data . 


The  2-hr  averaging  time  is  chosen  because  no  significant 
spectral  peaks  are  found  to  occur  for  periods^  less  than  10  hr. 

A  2-hr  averaging  time  guarantees  that  the  computed  PSD  is 
meaningful  for  periods  greater  than  4  hr,  i.e.,  the  Nyquist 
sampling  period.  Linear  interpolation  is  selected  for  data 
synthesis,  based  on  a  comparison  with  two  alternative  techniques: 


*Non-stationarity  occurs  when  the  length  of  the  time  series 
is  short  compared  to  the  period  of  the  low-frequency  components 
(Ref.  12). 

tThe  period  of  a  spectral  component  is  the  inverse  of  the  cyclical 
frequency . 
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placing  zeroes  in  the  gaps  and  ignoring  the  gaps  (i.e., 
"squeezing”  the  time  series).  Interpolation  preserves  the 
low-frequency  content  of  the  data,  and  does  not  introduce  dis¬ 
continuities  into  the  time  series  or  associated  spurious  fre¬ 
quencies  into  the  spectrum.  Detrending  is  an  important  aspect 
of  the  PSD-computation  procedure  in  cases  where  the  TD  values 
in  May  1979  and  May  1980  differ  significantly.  If  the  TD  trend 
is  not  removed,  the  PSD  and  ACF  are  distorted  by  the  occurrence 
of  high  energy  near  zero  frequency. 

ACFs  and  PSDs  computed  by  the  aforementioned  procedure 
are  presented  in  Figs.  3.3-1  to  3.3-3.  The  time  series  all 
contain  spectral  peaks  at  the  24-hr  and  12-hr  periods  (0.042 
cycle/hr  and  0.083  cycle/hr,  respectively,  in  the  PSD  plots). 
The  24-hr  diurnal  peak  is  10  to  20  dB  below  the  zero- frequency 
spectral  level  and,  in  some  cases,  is  more  than  10  dB  above 
the  spectral  level  at  adjacent  frequencies  (e.g.,  see  Fig. 
3.3-3b).  The  ACF  generally  decreases  with  increasing  time 
lag,  and  contains  a  sinusoid-like  component  with  a  period  of 
24  hr. 

The  width  of  the  24-hr  "hump"  in  the  PSD  depends  on 
the  seasonal  variation  (i.e.,  modulation)  of  the  diurnal-cycle 
amplitude.  The  hump  would  be  a  "spike"  if  the  diurnal  cycle 
were  a  pure  sinusoid  with  constant  amplitude.  The  12-hr  com¬ 
ponent  in  the  PSD  may  be  partly  due  to  the  bimodal  diurnal 
cycle  identified  in  Section  3.2.3.  However,  such  a  component 
is  expected  to  arise  even  for  a  uniraodal  diurnal  cycle,  if  the 
cycle  is  not  a  pure  amplitude-modulated  sinusoid.  Indeed,  the 
12-hr  component  is  the  first  harmonic  of  the  fundamental  24-hr 
component . 

The  slope  of  the  ACF  differs  significantly  for  the 
various  time  series.  The  slope  is  governed  largely  by  the 
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Figure  3.3-3  Autocorrelation  Functions  and  Power 
Spectral  Densities  for  TDZ 
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nature  of  the  low-frequency  (seasonal)  variations.  Because 
the  seasonal  variations  may  differ  from  year  to  year  and  the 
ACF  is  based  on  one  year  of  data,  the  ACFs  are  not  particularly 
useful  for  comparison  purposes.  In  other  words,  the  time  series 
are  not  strictly  stationary.  Some  improvement  in  stationarity 
can  be  achieved  by  high-pass  filtering  the  data  to  remove  seasonal 
variations  which  can  not  be  removed  by  simple  detrending.  A 
comparison  of  ACFs,  with  and  without  filtering,  is  presented 
in  Fig.  3.3-4  to  illustrate  this  point.  The  ACF  with  filtering 
crosses  zero  at  a  time  lag  of  approximately  3  days,  which  is 
consistent  with  U.S.  Coast  Guard  analyses  of  one-month  segments 
of  low-density  data  (e.g.,  Ref.  16).  Filtering  is  not  performed 
in  the  ACF/PSD  analyses  in  this  study,  to  permit  a  comparison 
between  the  actual  spectral  levels  at  zero  and  diurnal  frequencies. 
However,  it  should  be  kept  in  mind  that  the  relative  spectral 
levels  may  be  different  in  other  years. 


R63B15 


a)  No  Filter  (Same  as 
Fig.  3 . 3-2b) 


b)  High-Pass  Filter  (Cut-Off 
Period  =  30  days) 


Figure  3.3-4  Effect  of  High-Pass  Filtering  on  the  ACF 
(TDY  for  Dunbar/LC-204/2 ) 
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3.3.2  Correlation  of  TD  Pairs 


Relationships  among  the  various  TDs  are  identified  by 
correlating  pairs  of  time  series.  The  3-day  smoothed  seasonal 
time  series  presented  in  Section  3.2.1  are  used  for  this  purpose. 

Correlation  coefficients  (p)  for  TD  pairs  arranged  by 
site  are  presented  in  Tables  3.3-1  to  3.3-3.  Of  particular 
interest  is  the  correlation  of  redundant  TDs,  i.e.,  those 
measured  at  the  same  site  with  different  receivers.  Redundant 
TDY  data  collected  with  the  two  LC-204  receivers  at  DeTour 
exhibit  a  correlation  coefficient  of  0.94  and,  therefore,  are 
nearly  identical.  However,  the  redundant  TDY  data  collected 
with  the  BRN-5  and  LC-204  receivers  at  Dunbar  are  less  highly 
correlated  (p  =  0.80).  The  smaller  correlation  may  be  caused  by 
differences  in  the  Magnavox  and  Internav  signal-processing 
techniques  or  notch-filter  settings.  The  low  correlation  be¬ 
tween  the  TDZ  data  collected  with  the  two  LC-204  receivers  at 
Iroquois  (p  =  0.61)  suggests  that  at  least  one  of  these  receivers 
may  not  be  functioning  properly. 

Correlation  coefficients  for  TD  pairs  arranged  by 
station  are  presented  in  Tables  3.3-4  to  3.3-6.  These  results 
reveal  that  TDY  data  for  the  three  sites  are  highly  correlated 
(p  =  0.87  to  0.94).  This  high  correlation,  together  with  the 
low  correlation  for  TDX  and  TDZ,  suggests  that  TDY  variations 
are  caused  by  a  common  non-propagation  mechanism.  The  mechanism 
is  necessarily  associated  with  either  the  TDY  channel  of  the 
Austron-5000  receiver  at  the  SAM  or  with  transmitter  Y  itself. 

The  relationship  between  TDY  and  temperature  variations  is 
explored  in  Section  3.3.3. 


TABLE  3.3-1 

CORRELATION  COEFFICIENTS  FOR  DETOUR 


LC-204/1 

LC-204/2 

TDX  TDY 

TDY  TDZ 

LC-204/1 

TDX 

TDY 

1 '  "v>. 

l  l 

l  l 

LC-204/2 

TDY 

TD2 

+  1.00 

+0.05  +1.00 

TABLE  3.3-2 

CORRELATION  COEFFICIENTS  FOR  DUNBAR 


BRN-5/1 

LC-204/2 

H9 

TDX  -  TDY 

TDY  TDZ 

BRN-5/1 

i 

+  1.00 

+0.89  +1.00 

1  ! 

1  1 

LC-204/2 

■ 

Gnu 

9  1 

TABLE  3.3-3 

CORRELATION  COEFFICIENTS  FOR  IROQUOIS 


LC-204/1 

LC-204/2 

TDX  TDZ 

TDY  TDZ 

LC-204/1 

TDX 

TDZ 

+  1.00 

+0.41  +1.00 

l  l 

i  i 

LC-204/2 

+  1.00 

+0.71  +1.00 
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TABLE  3.3-4 

CORRELATION  COEFFICIENTS  FOR  TDX 


DETOUR 

DUNBAR 

IROQUOIS 

LC-204/1 

BRN-5/1 

LC-204/1 

DETOUR 

LC-204/1 

+  1.00 

- 

- 

DUNBAR 

BRN-5/1 

+0.22 

+  1.00 

- 

IROQUOIS 

LC-204/1 

+0.18 

+0.56 

+  1.00 

TABLE  3.3-5 

CORRELATION  COEFFICIENTS  FOR  TDY 


DETOUR 

DUNBAR 

IROQUOIS 

LC- 2t)4/2 

LC-204/2 

LC-204/2 

DETOUR 

DUNBAR 

IROQUOIS 

LC-204/2 

LC-204/2 

LC-204/2 

+1.00 

+0.94 

+0.87 

+  1.00 

+0.88 

+  1.00 

TABLE  3.3-6 

CORRELATION  COEFFICIENTS  FOR  TDZ 


DETOUR 

DUNBAR 

IROQUOIS 

LC-204/2 

LC-204/2 

LC-204/2 

DETOUR 

LC-204/2 

+  1.00 

- 

- 

DUNBAR 

LC-204/2 

+0.10 

+  1.00 

- 

IROQUOIS 

LC-204/2 

+0.07 

-0.17 

+  1.00 

Each  correlation  coefficient  presented  above  is  a 
measure  of  the  relationship  between  two  seasonal  time  series. 

It  is  also  of  interest  to  determine  the  short-term  correlation, 
by  employing  a  more  refined  time  scale.  The  coherency  spectrum , 
which  in  effect  shows  the  correlation  coefficient  as  a  function 
of  frequency,  is  very  useful  in  this  regard.  The  coherency 
spectrum  is  related  to  the  cross  PSD,  which  in  turn  is  computed 
as  the  Fast  Fourier  Transform  of  the  cross-correlation  function 
(Ref.  12).  Examples  of  the  cross-correlation  function  and 
coherency  spectrum  are  given  in  Fig.  3.3-5  for  TDX  vs  TDY  data 
from  the  Dunbar/BRN-5/1  receiver.  The  coherency  at  zero  fre¬ 
quency  is  approximately  0.89,  the  same  as  the  correlation  co¬ 
efficient  previously  given  in  Table  3.3-2.  However,  the  co¬ 
herency  decreases  rapidly  and  is  less  than  0.50  for  frequencies 
greater  than  0.01  cycle/hr  (i.e.,  periods  less  than  9  days). 

Based  on  these  results,  TD/TD  relationships  are  expected  to  be 
more  pronounced  on  a  seasonal  rather  than  diurnal  time  scale. 

The  value  of  the  cross-correlation  function  at  zero  time  lag 
(0.75;  see  Fig.  3.3-5)  can  be  interpreted  as  an  "ensemble  average" 
correlation  coefficient  over  all  frequencies.  However,  it 
should  be  recognized  that  this  average  value  is  weighted  by 
the  energy  at  each  frequency  and,  thus,  is  closer  to  the  zero- 
frequency  coherency  than  otherwise  expected. 


3.3.3  Correlation  of  TPs  with  Temperature 
and  Refractive  Index 


Meteorological  data  from  the  National  Weather  Service 
Station  at  Sault  Sainte  Marie,  Michigan  were  provided  to  TASC 
on  magnetic  tape  by  the  U.S.  Coast  Guard.  The  data  include 
the  dry-bulb,  wet-bulb,  and  dew-point  temperatures  and  the 
relative  humidity,  all  recorded  at  3-hr  intervals.  The  3-day 
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Figure  3.3-5  Cross-Correlation  Function  and  Coherency 
Spectrum  for  TDX  vs  TDY  Data  (Dunbar/ 

BRN-5/1 ) 

smoothed  dry-bulb  temperature  is  shown  in  Fig.  3.3-6  and  refrae- 
tivity",  computed  by  the  formula  in  Refs.  3  and  8,  is  presented 
in  Fig.  3.3-7. 

Correlation  coefficients  for  TDs  vs  temperature  and 
TDs  vs  refractivity  are  given  in  Table  3.3-7.  These  results 
indicate  a  relatively  strong  negative  correlation  (p  =  -0.77 
to  -0.85)  between  TDY  and  temperature  for  all  sites  and  receivers. 
Because  the  TD/temperature  correlation  is  generally  weaker  for 
TDY  than  for  TDX  and  TDZ,  it  appears  that  the  temperature  effect 
is  not  related  to  chain-wide  propagation  conditions.  This 
contention  is  further  justified  by  the  TD/refractivity  correla¬ 
tion  coefficients  presented  in  Table  3.3-7,  which  are  all 
less  than  or  equal  to  0.61  in  absolute  value.  The  order  of 
the  TD/refractivity  correlation  coefficients  is  the  same  as 
for  the  TD/temperature  correlation  coefficients,  because  of 
the  dependence  of  refractivity  on  temperature.  However,  the 
fundamental  parameter  is  clearly  temperature. 

-  ^ 

*Re f ractivi ty  =  (refractive  index  -  1.0)  *  10  . 
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Figure  3.3-7 


Refractivity  Seasonal  Time  Series 


J 


TABLE  3.3-7 

CORRELATION  OF  TDs  WITH  TEMPERATURE  AND  REFRACTIVITY 


SITE 

RECEIVER 

TYPE/ 

NUMBER 

TIME 

DIFFERENCE 

CORRELATION  COEFFICIENTS 

TD  vs 

TEMPERATURE 

TD  vs 

REFRACTIVITY 

DeTour 

LC-204/1 

X 

-0.22 

-0.14 

Y 

-0.77 

-0.55 

DeTour 

LC-204/2 

Y 

-0.84 

-0.57 

Z 

-0.19 

-0.24 

Dunbar 

BRN-5/1 

X 

-0.70 

Y 

-0.80 

-0.61 

Dunbar 

LC-204/2 

Y 

-0.85 

-0.61 

Z 

+0.54 

+0.35 

Iroquois 

LC-204/1 

X 

-0.76 

-0.56 

Z 

-0.51. 

-0.31  | 

Iroquois 

LC-204/2 

Y 

-0.80 

-0.60 

Z 

_ 

-0.66 

-0.51 

The  correlation  between  TDY  and  temperature  is  further 
explored  by  computing  the  cross-correlation  function  and  co¬ 
herency  spectrum  in  the  case  of  the  I roquois/LC-204/2  data 
(see  Fig.  3.3-8).  The  coherency  is  approximately  0.80  at  zero 
frequency,  consistent  with  the  correlation  coefficient  pre¬ 
sented  in  Table  3.3-7.  However,  the  coherency  falls  off  rapidly 
with  increasing  frequency,  suggesting  that  the  effect  of  tempera¬ 
ture  is  primarily  seasonal.  The  U.S.  Coast  Guard  has  computed 
TD/temperature  correlation  coefficients  based  on  14-day  smoothing 
(Ref.  21).  The  results  show  larger  TDZ/temperature  correlation 
coefficients  for  the  Dunbar  and  Iroquois  sites  than  those  com¬ 
puted  based  on  3-day  smoothing  in  the  present  study.  The  dis¬ 
crepancy  is  consistent  with  the  coherency  spectrum. 
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Figure  3.3-8  Cross-Correlation  Function  and 
Coherency  Spectrum  for  TDY 
( Iroquois/LC-204/2 )  vs  Temperature 
Data 

3. 4  LPA  STATISTICS 

The  LPA  data  contain  information  regarding  the  propaga 
tion  variations  which  occur  along  transmi tter- to-SAM  signal 
paths.  However,  extraction  of  the  information  requires  careful 
removal  of  clock  phase  drift  caused  by  differences  in  the  fre¬ 
quencies  of  the  secondary  and  master  cesium  standards. 

To  illustrate  this  point,  the  daily  sum  of  the  LPAs 
is  presented  in  Fig.  3.4-1  for  each  secondary  transmitter. 
Inspection  of  the  results  for  transmitter  X  might  lead  to  the 
conclusion  that  propagation  variations  are  larger  in  September 
than  in  January.  However,  a  graph  of  the  cumulative  LPA  sum 
vs  time  (see  Fig.  3.4-2)  shows  that  the  large  difference  in 
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daily  LPA  sums  is  caused  by  a  difference  in  clock  drift  rates 
for  the  two  months.  The  drift  rate  in  September  is  approxi¬ 
mately  0.17  psec/day,  but  after  clock  adjustments  made  in 
November  the  drift  rate  is  only  -0.02  psee/day.  Therefore, 
the  daily  LPA  sums  for  September  shown  in  Fig.  3.4-la  are 
primarily  associated  with  clock  synchronization ,  rather  than 
propagation.  Separation  of  synchronization  and  propagation 
effects  is  beyond  the  scope  of  the  current  effort. 


3.5  COMPARISON  OF  THEORETICAL  AND 
OBSERVED  I'D  VARIATIONS 

Signal  propagation  theory  is  employed  in  Ref.  3  to 
estimate  the  magnitude  of  St.  Marys  River  Loran-C  TD  grid  in¬ 
stability.  Because  the  mini-chain  coverage  area  is  less  than 
60  nm  in  diameter,  the  propagation  medium  is  assumed  to  be 
homogeneous.  This  simpli fying . assumption  is  consistent  with 
the  goal  of  estimating  temporal  TD  variations,  although  it  is 
not  necessarily  appropriate  for  a  model  of  spatial  variations 
(Refs.  1  and  2). 

The  sensitivity  of  TDs  to  temporal  variations  in 
three  propagation  parameters  --  atmospheric  refractive  index, 
vertical  lapse  rate  of  refractive  index,  and  ground  conduc¬ 
tivity  --  is  analyzed  in  Ref.  3.  It  is  shown  that  TD  sensi¬ 
tivity  to  refractive  index  and  vertical  lapse  rate  variations 
is  proportional  to  the  "Double  Range  Difference" 

DRD  =  <RS  -  Rm)  -  <R;  -  R;>  (3.5-1) 

where  Rg  and  Rffi  are  the  secondary-to-site  master-to-site  ranges, 
respectively,  and  R^  and  R^  are  the  corresponding  ranges  to 
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the  SAM.  DRD  contours  are  hyperbolas,  as  shown  in  Figs.  3.5-1 
to  3.5-3.  The  contour  of  zero  sensitivity  (DRD  =  0)  passes 
through  the  SAM  location,  because  the  SAM  maintains  its  TD  at 
a  constant  level.  TD  sensitivity  to  conductivity  variations 
is  given  by  a  Modified  Double  Range  Difference 

Modified  DRD  =  (R  a  -  Rma)  -  (R'a  -  R’a)  (3.5-2) 

w}  Uk  3  111 

where  the  exponent  "a"  approximately  equals  0.5.  This  modi¬ 
fication  accounts  for  the  nonlinearity  in  the  secondary  phase 
delay  vs.  range  function  for  ranges  less  than  100  nm  (Ref.  15). 


Figure  3.5-1 


TDX  Temporal  Variation  Sensitivity 
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Figure  3.5-2  TDY  Temporal  Variation  Sensitivity 


The  theoretical  TD  variations  associated  with  particu 
variations  in  the  propagation  parameters  are  indicated  below: 

JL 

•  Refractive  Index  --  6x10  psec/nm'  per 
N-uni t 

•  Vertical  Lapse  Rate  --  3x10*^  psec/nm  per 
0.01  a-units 

_  3 

•  Conductivity  --  4x10  psec/nm  per  doubling 
of  conductivity. 


*TD  variation  per  nm  of  DRD.  For  conductivity,  a  linear 
approximation  is  used  to  permit  the  DRD  to  replace  the 
Modified  DRD. 
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Figure  3.5-3  TDZ  Temporal  Variation  Sensitivity 


The  maximum  DRD  for  the  St.  Marys  River  Loran-C  data  collec¬ 
tion  sites  is  35  nm  (see  Figs.  3.5-1  to  3.5-3).  Therefore, 
the  maximum  seasonal  TD  variation  expected  for  a  50  N-unit 
refractive  index  variation  (see  Fig.  3.3-7)  is  only  0.010  psec. 
Similarly,  the  maximum  seasonal  TD  variation  expected  for  a 
0.25  a-unit  vertical  lapse  rate  variation  (Refs.  3  and  9)  is 
0.026  psec.  Clearly,  the  ±0.2-(Jsec  seasonal  TD  variations 
shown  in  Section  3.2.1,  if  propagation-related,  are  necessarily 
caused  by  conductivity  variations.  The  maximum  TD  variation 
expected  for  a  doubling  of  conductivity  from  winter  to  summer 
(Ref.  3)  is  0.14  psec  for  the  three  data  collection  sites. 
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Although  the  maximum  magnitude  of  the  observed  TD 
variations  is  consistent  with  propagation  theory,  the  relative 
variations  of  the  individual  time  series  are  not  consistent 
with  the  DRD.  This  is  illustrated  in  Fig.  3.5-4,  a  comparison 
of  the  seasonal  time  series  and  the  corresponding  DRDs .  The 
following  are  among  the  obvious  discrepancies: 

•  Dunbar  TDY  variations  are  0.189  psec 
peak- to-peak ,  but  the  DRD  equals  zero 

•  DeTour  TDZ  variations  are  only  0.061 
psec  peak- to-peak ,  but  the  DRD  equals 
+35  nm 

•  DeTour  and  Iroquois  TDY  variations  are 
in  the  same  direction,  but  the  DRDs  are 
of  opposite  sign  (-25  nm  vs  +8  nm) 

•  DeTour  and  Iroquois  TDX  variations  differ 
significantly,  but  the  DRDs  are  both 
equal  to  +11  nm. 

For  the  data  collection  sites  considered,  the  Modified  DRD  is 
nearly  proportional  to  the  DRD  and,  therefore,  does  not  contri 
bute  to  resolving  the  discrepancies. 

It  is  possible  that  seasonal  conductivity  variations 
are  spatially  heterogeneous,  thereby  violating  the  assumptions 
inherent  in  the  DRD  model.  However,  it  is  more  likely  that 
the  observed  TD  variations  are  caused  by  a  transmitter-  and/or 
receiver-related  mechanism,  which  dominates  or  disguises  the 
propagation  mechanism. 
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Figure  3.5-4  Comparison  of  Seasonal  TD  Time  Series 


The  U.S.  Coast  Guard  has  conducted  low-density  data 
analyses,  which  closely  parallel  the  high-density  analyses 
presented  in  this  report.  The  low-density  approach  is  motivated 
by  the  desire  of  the  U.S.  Coast  Guard  to  have  a  quick- response , 
in-house  analysis  capability.  Specifically,  the  reduced  data 
volume  has  enabled  monthly  internal  reports  to  be  issued,  which 
summarize  the  analysis  results  (e.g..  Ref.  6).  A  statistical 
evaluation  of  the  low-density  approach  is  presented  in  this 
chapter  to  validate  the  monthly  reports  and  to  establish  an 
efficient  Loran-C  monitoring  strategy  for  harbors. 


4.2  COMPARISON  OF  EDITED  LOW-DENSITY  AND 
HIGH-DENSITY  DATA 

Potential  errors  in  the  low-density  TD  data,  compared 
to  the  high-density  data,  can  be  categorized  as  follows: 

•  Relatively  small  errors  associated  with 
computing  mean  TD  values  based  on  a  reduced 
sample  size 

•  Relatively  large  errors  resulting  from 
the  inability  to  detect  outliers  in  short 
segments  of  a  time  series. 


The  first  error  category  is  addressed  in  this  section,  by 
employing  low-density  data  extracted  from  the  edited  high- 
density  data  base.  Such  low-density  data  would  result  if 


4-1 


there  were  no  outliers  in  the  measured  time  series.  The  con¬ 
tribution  of  outliers  is  explored  in  Section  4.3. 

"Generic"  low-density  sampling  strategies  are  employed 
to  investigate  the  effect  of  sample  size.  Each  generic  strategy 
is  defined  by  three  parameters: 

•  Number  of  equally-spaced  time  windows 
per  day,  during  which  15-min  samples  are 
recorded 

•  Width  of  time  windows 

•  Phase  of  time  windows,  i.e.,  the  time-of- 
day  at  the  beginning  of  the  first  window. 

The  nominal  U.S.  Coast  Guard  sampling  strategy  involves  two 
time  windows  of  width  equal  to  1  hr,  beginning  at  0  hr  and  7  hr 
GMT.  Because  the  windows  are  not  equal ly- spaced  (i.e.,  12  hr 
apart),  this  strategy  is  addressed  separately  from  the  generic 
strategies . 

A. 2.1  Seasonal  Variations 

The  utility  of  the  low-density  data  in  quantifying 
seasonal  TD  variations  depends  on  the  difference  between  the 
low-  and  high-density  3-day  smoothed  time  series.  The  rms 
differences  (i.e,  errors)  for  the  generic  sampling  strategies 
are  presented  as  a  function  of  window  width  in  Fig.  4.2-la  and 
as  a  function  of  window  phase  in  Fig.  4.2-lb.  These  results 
are  obtained  by  processing  TDY  data  from  the  lroquois/LC-204/2 
receiver,  which  is  selected  to  represent  worst-case  seasonal 
variations.  A  direct  comparison  of  the  high-density  and  U.S. 
Coast  Guard  low-density  time  series  is  given  in  Fig.  4.2-lc. 
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Figure  4.2-1 


Accuracy  of  Edited  Low-Density  Seasonal  Time 
Series  (TDY  for  1 roquoi s/LC-204/2 ) 
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The  following  conclusions  can  be  drawn,  regarding  the 
accuracy  of  the  low-density  seasonal  time  series: 

•  TD  errors  are  less  than  0.01  psec  rms 
for  all  strategies  considered 

•  Four  1-hr  windows  are  more  effective 
than  two  2-hr  windows,  although  they  con¬ 
tain  the  same  number  of  samples 

•  TD  errors  are  more  sensitive  to  window 
phase  for  two  windows  than  for  four  windows 

•  The  U.S.  Coast  Guard  sampling  strategy 
could  be  improved  only  slightly  (0.0081 
psec  error  to  0.0067  psec  error)  by  plac¬ 
ing  the  windows  at  0  hr  and  12  hr  GMT. 

If  TD  errors  are  less  than  0.01  psec  rms,  as  is  the  case  for 
the  U.S.  Coast  Guard  sampling  strategy,  seasonal  variations 
can  be  adequately  identified.  Although  a  significant  reduc¬ 
tion  in  TD  errors  can  be  achieved  by  using  four  1-hr  windows 
instead  of  two,  the  additional  detail  does  not  contribute 
significantly  to  the  quantification  of  seasonal  variations. 

4.2.2  Monthly  Variations 

The  utility  of  the  low-density  data  in  quantifying 
monthly  TD  variations  depends  on  the  difference  between  the 
low-  and  high-density  6-hr  smoothed  time  series.  The  rms  TD 
errors  for  the  j-eneric  sampling  strategies  are  presented  as  a 
function  of  window  phase  for  July  and  January,  in  Figs.  4.2-2a 
and  4.2-2b,  respectively.  These  graphs  are  based  on  TDX  data 
from  the  Iroquois/LC-204/1  receiver,  which  exhibits  worst-case 
TD  variations  for  January.  A  comparison  of  the  high-density 
and  U.S.  Coast  Guard  low-density  time  series  for  January  is 
presented  in  Fig.  4.2-2c. 
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F igure  4.2-2 


The  low-density  monthly  time  series  is  somewhat  less 
accurate  than  the  low-density  seasonal  time  series  but,  never¬ 
theless,  is  in  error  by  less  than  0.012  psec  rms .  Consequently, 
the  U.S.  Coast  Guard  sampling  strategy  is  judged  to  be  adequate 
for  measuring  monthly  TD  variations.  The  large  sensitivity  of 
TD  errors  to  window  phase,  especially  for  the  January  data 
(see  Fig.  4.2-2b),  is  apparently  not  related  to  the  diurnal  TD 
cycle.  Indeed,  the  TDX  diurnal  cycle  for  January  for  the 
lroquois/LC-204/1  receiver  is  less  than  0.005  psec  peak-to- 
peak  (see  Fig.  3.2-l6c).  Rather,  the  two  minima  in  the  TD- 
error  vs  window-phase  graph  arise  from  the  assumed  definition 
of  a  monthly  time  series.  Specifically,  the  best  estimate  of 
the  6-hr  means  is  obtained  by  sampling  in  the  center  of  the 
6-hr  intervals,  i.e.,  at  3  hr,  9  hr,  15  hr,  and  21  hr  GMT. 

4.2.3  Diurnal  Variations 

Intuitively,  it  is  difficult  to  accurately  measure 
the  peak-to-peak  value  of  the  diurnal  TD  cycle  by  using  only 
two  1-hr  windows.  The  phase  of  the  windows  is  a  critical 
parameter  in  this  case,  more  so  than  in  the  measurement  of 
seasonal  and  monthly  TD  variations. 

The  error  in  the  peak-to-peak  value  of  the  low-density 
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diurnal  cycle  is  presented  as  a  function  of  window  phase  for 
July  and  January,  in  Figs.  4.2-3a  and  4.2-3b,  respectively. 

TDZ  data  from  the  lroquois/LC-204/2  receiver  are  selected  for 
the  analyses,  because  they  exhibit  a  distinct  diurnal  cycle. 

The  high-density  and  U.S.  Coast  Guard  low-density  diurnal  cycles 
are  compared  for  July  and  January,  in  Figs.  4.2-3c  and  4.2-3d, 
respectively . 


*For  two  1-hr  windows,  the  "diurnal  cycle”  consists  of  two 
points,  each  computed  as  described  in  Section  3.2.3. 
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The  following  conclusions  are  drawn  regarding  the 
accuracy  of  the  low-density  diurnal  cycle: 


•  TD  cycle  errors  are  more  sensitive  to 
window  phase  in  January  than  in  July, 
because  the  diurnal  cycle  amplitude  is 
larger  in  January 

•  The  optimal  window  phase  is  different  in 
July  than  in  January,  because  the  phase 
of  the  diurnal  cycle  is  different  in  the 
two  months 

•  TD  cycle  errors  for  January,  range  from 
less  than  0.001  psec  to  0.020  psec,  de¬ 
pending  on  the  choice  of  window  phase 

•  Sensitivity  to  window  phase  is  greatly 
reduced  by  using  four  windows  instead  of 
two 

•  The  "peak"  in  the  July  diurnal  cycle  and 
the  "trough"  in  the  January  diurnal  cycle 
are  not  observable  with  the  U.S.  Coast 
Guard  sampling  strategy,  thereby  result¬ 
ing  in  a  50-percent  error  in  the  measured 
cycle  amplitudes. 


It  is  important  to  recognize  that  the  optimal  window  phase  is 
generally  different  for  the  other  TD  time  series  (see  Figs. 
3.2-13  to  3.2-18). 

Clearly,  the  U.S.  Coast  Guard  low-density  data  should 
not  be  employed  for  diurnal  analyses.  Use  of  the  data  to  iden¬ 
tify  spectral  components  for  periods  less  than  or  equal  to 
24  hr  leads  to  an  aliasing  error  (Ref.  13).  To  minimize  this 
problem,  it  is  recommended  that  four  or  more  time  windows  per 
day  be  employed  in  the  Loran-C  harbor  monitor  program. 
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4.3  EFFECT  OF  OUTLIERS  ON  LOW-DENSITY  DATA 

The  high-density  editing  procedure  described  in  Chap¬ 
ter  2  consists  of  three  steps:  signal-quality  check  (Editing 
Mark  2),  tolerance  check  (Editing  Mark  3),  and  outlier  editing 
(Editing  Mark  6).  A  signal-quality  check  is  performed  in  the 
microprocessors  of  the  low-density  data-recording  system,  and 
a  tolerance  check  is  performed  by  the  U.S.  Coast  Guard  upon 
receipt  of  the  data.  However,  the  ability  to  detect  outliers 
is  limited  by  the  short  (1-hr)  contiguous  time  periods  during 
which  15-min  samples  are  available.  In  particular,  there  are 
not  enough  samples  to  compute  a  reliable  daily  trend  line  (see 
Section  2.2.6). 

Outliers  occur  in  a  small  percentage  of  the  1-hr  win¬ 
dows:  only  1.9  to  10.8  percent  of  the  high  density-data  are 
outliers,  and  only  one-twelfth  (8.3  percent)  of  the  high-den¬ 
sity  data  are  involved  in  the  low-density  analyses.  Therefore, 
a  comparison  of  time  series,  analogous  to  the  comparisons  pre¬ 
sented  in  Section  4.2,  is  not  particularly  useful  for  showing 
the  effect  of  outliers.  The  following  alternative  approach  is 
selected  and  applied  to  the  TDY  time  series  for  the  DeTour/LC-204/1 
receiver : 

•  Partition  the  time  series  into  1-hr  inter¬ 
vals  and  consider  only  those  intervals 
which  contain  outliers 

•  Calculate  the  differences  (errors)  between 
the  1-hr  TD  means  computed  with  and  without 
outliers 

•  Construct  a  histogram  of  the  TD  errors. 
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The  histogram  shows  the  distribution  of  all  potential  errors 
encountered  in  a  low-density  monthly  time  series,  due  to  the 
presence  of  outliers.  However,  very  few  of  the  errors  are 
actually  encountered,  because  very  few  1-hr  intervals  with 
outliers  are  involved  in  the  low-density  data. 

Histograms  are  presented  in  Fig.  4.3-1  for  three  defini¬ 
tions  of  outliers:  Editing  Mark  6  =  "1",  "1"  to  "4",  or  "1" 
to  "6".  It  is  unlikely  that  outliers  with  Editing  Mark  6  = 

"1"  can  be  detected  in  the  low-density  data.  However,  the 
potential  TD  errors  introduced  by  these  outliers  are  typically 
less  than  0.04  psec,  not  a  major  concern.  Outliers  with  Edit¬ 
ing  Mark  6  =  "2"  to  "4"  can  be  detected  by  a  careful  inspec¬ 
tion  of  the  low-density  data,  except  when  multiple  outliers 
occur  in  a  1-hr  window.  Such  an  inspection  is  warranted  be¬ 
cause  the  TD  errors  can  be  as  large  as  0.10  psec.  As  indicated 
in  Fig.  4.3-lc,  outliers  with  Editing  Mark  6  =  "5”  or  "6"  can 
cause  TD  errors  as  large  as  0.2  psec.  Fortunately,  these  gross 
outliers  are  easily  detected  and  removed  from  the  low-density 
data.  The  implication  of  the  skewed  histogram  in  Fig.  4.3-lc 
is  that  gross  outliers  consistently  lie  below  the  daily  trend 
line,  suggesting  a  possible  anomaly  in  the  DeTour/  LC-204/1 
receiver . 


4.4  ASSESSMENT  OF  LOW- DENSITY  DATA  ANALYSIS  TECHNIQUES 

The  following  are  among  the  low-density  data  analysis 
techniques  used  by  the  U.S.  Coast  Guard  (Refs.  17  and  18): 


•  K-Factor  Analysis  --  used  to  predict 
relative  TD  variations  in  an  assumed 
homogeneous  chain  coverage  area 


*The  low-density  monthly  time  series  necessarily  consists 
of  1-hr  means,  rather  than  6-hr  means. 
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Figure  4.3-1  Effect  of  Outliers  on  Low-Density  Monthly  Time 
Series  (TDY  for  DeTour/LC-204/1 ) 
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•  Uniform  Propagation  Velocity  Model  -- 
used  to  extract  propagation  velocity  and 
common  SAM- induced  errors  from  the  TD 
data 

•  Radial  Position  Error  Estimation  --used 
to  project  propagation  velocity  varia¬ 
tions,  SAM-induced  errors,  and  measure¬ 
ment  noise  into  Loran-C  user  position 
errors 

•  TD/Temperature  Correlation  Analysis  -- 
used  to  confirm  the  TD/vertical  lapse 
rate  relationship  expected  by  the  U.S, 

Coast  Guard. 

TASC  believes  these  analysis  techniques  are  appropriately  and 
correctly  applied,  except  as  noted  below. 

K-factor  analysis  is  equivalent  to  the  Double  Range 
Difference  analysis  presented  in  Section  3.5.  The  technique 
is  based  on  the  assumption  that  the  propagation  delay  increases 
linearly  with  range.  If  conductivity  variations  are  more  impor¬ 
tant  than  vertical  lapse  rate  variations,  however,  the  depen¬ 
dence  is  nonlinear  in  range.  In  this  case,  a  modified  K-factor, 
analogous  to  the  Modified  Double  Range  Difference  (see  Eq . 
3.5-2),  is  a  more  accurate  index, 

A  similar  modification  to  the  geometry  matrix  of  the 
uniform  propagation  velocity  moocl  is  required  to  accommodate 
conductivity  variations.  It  is  also  cautioned  that  the  "common 
SAM-induced  error"  appearing  in  the  model  includes  a  zero-range 
bias  in  the  propagation  velocity  vs  range  function.  This  arti¬ 
ficial  bias  occurs  because  the  data  cover  a  limited  set  of 
ranges  and  the  assumed  velocity  vs  range  function  is  an  approxi¬ 
mation.  The  "bias"  varies  with  time,  in  conjunction  with  the 
modeled  propagation  velocity  variations. 
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It  is  shown  in  Ref.  1  that  the  cross-channel  position 
error  is  significantly  smaller  than  the  radial  (specifically, 
the  2d  rms)  position  error,  in  certain  areas  of  the  St.  Marys 
River  navigation  channel.  Therefore,  it  is  recommended  that 
both  the  cross-channel  and  radial  position  errors,  or  the  cross¬ 
channel  and  along-channel  position  errors,  be  included  in  com¬ 
parisons  with  the  channel  width  (Ref.  17). 

TD/temperature  correlation  analyses,  similar  to  those 
presented  in  Section  3.3.3,  are  an  important  part  of  the  U.S. 
Coast  Guard  low-density  approach.  It  is  recommended  that  the 
analyses  be  limited  to  correlation  coefficients.  The  cross¬ 
correlation  function  is  very  sensitive  to  the  nature  of  low- 
frequency  effects  and,  therefore,  cannnot  be  relied  on  for 
comparison  purposes.  Although  agreeing  that  TD/temperature 
relationships  must  be  identified,  TASC  is  not  in  agreement 
with  the  U.S.  Coast  Guard  contention  that  such  relationships 
are  primarily  a  vertical  lapse  rate  effect  (Ref.  17).  The 
theoretical  results  in  Refs.  3  and  9  show  clearly  that  the  TD 
variations  induced  by  vertical  lapse  rate  variations  are 
negligible  (<0.03  psec)  in  a  short-baseline  chain.  It  appears 
from  the  results  presented  in  this  report  that  the  high  TD/ 
temperature  correlation  observed  in  the  data  is  not  related  to 
propagation.  To  further  test  this  hypothesis,  it  is  recommended 
that  vertical  lapse  rate  data  from  the  National  Weather  Service 
Station  at  Sault  Sainte  Marie,  Michigan  be  correlated  against 
the  TD  data.  In  computing  the  vertical  lapse  rate  from  radio¬ 
sonde  data,  it  should  be  kept  in  mind  that  the  results  may  be 
sensitive  to  the  altitudes  selected  for  the  computation  (Ref. 
20).  Various  altitudes  in  addition  to  the  surface/l-km  stan¬ 
dards  should  be  tested. 


V 


CONCLUSIONS  AND  RECOMMENDATIONS 


Results  of  a  detailed  analysis  of  St.  Marys  River 
temporal  Loran-C  data,  collected  by  the  U.S.  Coast  Guard  from 
May  1979  to  May  1980,  are  documented  in  this  report.  The 
following  objectives  are  met  by  the  data  analysis  effort: 


•  Quantification  of  Time  Difference  (TD) 
grid  instability  in  the  St.  Marys  River 
Loran-C  mini-chain,  based  on  high-density 
(15-min)  data 

•  Establishment  of  an  efficient  Loran-C 
monitoring  strategy  for  harbors,  based 
on  an  evaluation  of  the  U.S.  Coast  Guard 

low-density  data  analysis  approach. 


Data-organizing  and  -editing  capabilities,  needed  to  support 
the  analyses,  are  provided  by  the  TASC  Loran-C  Data-Base  Manage¬ 
ment  Software  developed  in  a  previous  effort  (Ref.  5).  Conclu¬ 
sions  and  recommendations  related  to  the  two  study  objectives 
are  presented  below. 


5.1  QUANTIFICATION  OF  ST.  MARYS  RIVER  LORAN-C 
TD  GRID  INSTABILITY 


The  following  principal  conclusions  are  drawn  from 
analyses  of  the  high-density  TD  time  series  data: 


*Low-density  data  consist  of  five  15-min  samples  recorded 
twice  daily. 
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•  Seasonal  TD  variations  (i.e.,  variations 
in  3-day  means  over  the  data  collection 
period)  are  less  than  0.4  psec  peak-to- 
peak  at  the  three  data  collection  sites, 
consistent  with  previous  observations 
(Refs.  2  and  3) 

•  TD  variations  during  monthly  periods  are 
larger  in  the  winter  than  in  the  summer 
(e.g.,  see  Fig.  5.1-1),  consistent  with 
signal  propagation  theory  (Refs.  3  and  9) 

•  Diurnal  TD  cycles  are  generally  less 
than  0.02  psec  peak- to-peak ,  which  is 
small  relative  to  longer-term  variations 
(e.g. ,  see  Fig .  5.1-2) 

•  Power  Spectral  Densities  (PSDs)  computed 
from  the  data  show  that  the  diurnal  TD 
cycle  is  10  to  20  dB  below  the  zero- 
frequency  spectral  component  (e.g.,  see 
Fig.  5.1-3) 

•  TDZ  data  measured  at  the  Pt.  Iroquois 
site  differ  significantly  for  the  two 
Internav  LC-204  receivers  (the  correla¬ 
tion  coefficient,  p,  is  0.61) 

•  TDY  data  from  the  three  sites  are  highly 
correlated  (p  =  0.87  to  0.94) 

•  TDY  data  exhibit  a  relatively  high  nega¬ 
tive  correlation  with  Sault  Sainte  Marie 
temperature  data  (p  =  -0.77  to  -0.85). 


A  comparison  of  the  data  analysis  results  to  the  theoretical 
predictions  given  in  Ref.  3  shows  that  the  maximum  observed  TD 
variation  is  consistent  with  expected  seasonal  conductivity 
variations.  However,  the  relative  variations  of  the  individual 
time  series  are  not  consistent  with  propagation  theory  for  a 
homogeneous  medium. 


It  is  possible  that  seasonal  conductivity  variations 
are  spatially  heterogeneous,  thereby  violating  the  assumptions 
made  in  Ref.  3.  However,  it  appears  more  likely  that  the  ob¬ 
served  TD  variations  are  caused  by  a  transmitter-  and/or 
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receiver- related  mechanism,  which  disguises  the  propagation 
effect.  It  is  important  to  establish  the  nature  of  this 
mechanism,  because  of  its  potential  impact  on  the  Loran-C  har 
bor  monitor  program.  Therefore,  it  is  recommended  that  addi¬ 
tional  data  analyses  be  conducted,  to  search  for  possible 
cause/effect  relationships,  whether  they  be  propagation-  or 
transmitter/receiver- induced  (Ref.  19).  The  data  base 
described  herein  is  the  most  comprehensive  of  its  kind  ever 
assembled  and  provides  a  strong  nucleus  for  such  an  effort. 


5.2  EVALUATION  OF  LOW-DENSITY  DATA  ANALYSIS  APPROACH 


The  following  principal  conclusions  are  drawn  from 
the  evaluation  of  the  U.S.  Coast  Guard  low-density  data  analy 
sis  approach: 


•  Errors  in  the  low-density  seasonal  TD 
time  series,  relative  to  the  high-density 
time  series,  are  less  than  0.01  psec  rms 

•  A  reduction  in  the  low-density  seasonal 
time  series  errors  can  be  realized  by 
increasing  the  number  of  samples  collected 
per  day,  but  the  reduction  results  in 
only  a  slight  improvement  in  data  utility 

•  The  above  conclusions  also  hold  for  the 
low-density  monthly  time  series 

•  Significant  errors  (as  large  as  50  percent) 
result  if  the  low-density  data  are  used 

to  compute  the  peak-to-peak  value  of  the 
diurnal  TD  cycle 

•  The  low-density  d iurnal -eye le  errors  are 
very  sensitive  to  the  times-of-day  selected 
for  sampling 

•  The  optimal  sampling  times  for  diurnal -eye le 
resolution  generally  differ  for  different 
months  and  different  receiver/TD  combinations 


AD-AlOtS  07H 


UNCLASSIF 1ED 


ANALYTIC  SCIENCES  CORP  READING  -A 

*2T 0F  ST*  MAm  R,vc*  loran-c  time  OtrPERFNCE  GRID  — ETC  *  u> 
AUG  80  00T-CR-M-77-17B* 

WSCG-O-52-81  mi 


Moderate  data  outliers  are  not  a  problem 
in  the  low-density  data,  but  gross  outliers 
must  be  removed  by  inspection. 


Based  on  these  conclusions,  it  is  recommended  that  the  U.S. 

Coast  Guard  low-density  data  not  be  employed  for  diurnal  analyses. 
To  permit  diurnal  analyses  in  the  Loran-C  harbor  monitor  program, 
the  number  of  1-hr  time  windows  per  day  should  be  increased 
from  two,  to  at  least  four.  Furthermore,  it  must  be  recognized 
that  the  analysis  results  herein  pertain  specifically  to  the 
St.  Marys  River  mini-chain.  In  anticipation  that  different 
results  may  apply  to  a  long-baseline  chain,  it  is  recommended 
that  high-density  harbor  monitor  data  be  collected  during  both 
a  summer  month  and  a  winter  month  and  analyzed  prior  to  stan¬ 
dardizing  a  low-density  sampling  strategy. 

The  U.S.  Coast  Guard  low-density  data  analysis  techniques 
described  in  Ref.  17  are  judged  to  be  appropriately  and  correctly 
applied,  with  a  few  minor  exceptions  (see  Section  4. A).  However, 
TASC  is  not  in  agreement  with  the  U.S.  Coast  Guard  contention 
that  the  high  TD/temperature  correlation  observed  in  the  data 
is  primarily  a  vertical  lapse  rate  effect.  Based  on  signal 
propagation  theory,  the  vertical  lapse  rate  effect  is  expected 
to  be  negligible  in  a  mini-chain.  To  resolve  this  issue,  it 
is  recommended  that  vertical  lapse  rate  data  from  the  National 
Weather  Service  Station  at  Sault  Sainte  Marie,  Michigan  be 
correlated  against  the  TD  data. 


APPENDIX  A 

FORMATS  FOR  LORAN-C  DATA  SUPPLIED  TO  TASC 
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Formats  for  Loran-C  data  supplied  to  TASC  on  half-inch 
magnetic  tape  by  the  U.S.  Coast  Guard  are  presented  in  Tables 
A- 1  to  A-4.  The  record  lengths  are  80  characters  for  the  LC-204 
and  BRN-5  formats  and  40  characters  for  the  Austron-5000  format. 
All  characters  which  are  not  included  in  the  tables  are  "blanks." 
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TABLE  A- 1 

INTERNAV  LC-204  FORMAT  FOR  DATA  SUPPLIED  TO  TASC 


73484 


CHARACTERS 

DATA  ITEM 

■■■■■■■■■■■■  '1  fliTM***** 

1 

SWITCH 

Control  unit  switch  setting  (G  for 
valid,  B  for  invalid) 

2-4 

DAY 

Julian  day  (Greenwich  Mean  Time) 

5-6 

HR 

Hours 

7-8 

MIN 

Minutes 

9-10 

SEC 

Seconds 

11-14 

SNRECDM 

Master  Signal- to-Noise  Ratio  (SNR) 
or  Envelope-to-Cycle  Difference 
(ECD)  (See  char.  38) 

15-18 

SNRECD1 

Secondary  1  SNR  or  ECD 

20-28 

TD1 

Tine  difference  1  (0.1  nsec) 

33 

SITE 

Site  number 

34 

... 

(1) 

35-37 

DATAID 

Experiment  number 

38 

IDSNRECD 

Indicates  whether  SNR(0)  or  ECD(l) 
is  printed  this  record 

41-44 

SNRECD2 

Secondary  2  SNR  or  ECD 

45 

ECSM 

Master  Envelope  Command  Signal  (ECS) 
indicates  whether  no  cycle  jump  (0), 
positive  jump  (1),  or  negative  jump 
(2)  is  commanded  by  receiver 

46 

ECS1 

Secondary  1  ECS 

47 

ECS  2 

Secondary  2  ECS 

48 

SUFFS1G 

Sufficient  signal  indicator:  When 
converted  to  a  binary  number,  the  3 
digits  indicate  whether  master, 
secondary  1  and  secondary  2  signals 
are  insufficient  (0)  or  sufficient  ( 

50-58 


TD2 


Tine  difference  2  (0.1  nsec) 


TABLE  A- 2 

MAGNAVOX  AN/BRN-5  FORMAT  (SECONDARY)  FOR  DATA 
SUPPLIED  TO  TASC 


T3487 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

1 

1DENT 

Record  type  indicator:  secondary  1  (2) 
or  secondary  2  (3) 

2-4 

DAY 

Julian  day  (Greenwich  Mean  Time) 

5-6 

HR 

Hours 

7-8 

MIN 

Minutes 

9-10 

SEC 

Seconds 

12 

SITE 

Sice  number 

13-15 

DATA ID 

Experiment  number 

17 

SUFSIG 

Secondary  1  or  2  signal  status:  When  converted 
to  a  binary  number,  first  digit  indicates 
whether  (1)  or  not  (0)  signal  is  present  and 
second  digit  whether  (1)  or  not  (0)  station 
is  blinking 

19-27 

TD 

Time  difference  1  or  2  (0.1  nsec) 

29-37 

TOA 

Time  of  arrival  1  or  2  (0.1  nsec) 

39-42 

RMSTOA 

RMS  time  of  arrival  1  or  2  (nsec) 

4*4  -  4b 

ETE 

Secondary  1  or  2  envelope  time  estimate 
(0.1  usee).  Similar  to  ECD 

48-53 

SIGLEV 

Secondary  1  or  2  signal  level  (0.24414  mv) 

55-58 

FRQNT12 

Notch  filter  frequency  (0.1  KHz ) .  For 
notch  1  if  char.  1=2;  for  notch  2  if 
char.  1=3 

60-61 

RECATTEN 

Receiver  attenuation  (dB)  for  secondary 

1  or  2  signal 

63-66 

SNR 

Secondary  1  or  2  signal-to-noise  ratio  (0.1  dB) 

TABLE  A- 3 

MAGNAVOX  AN/BRN-5  FORMAT  (MASTER)  FOR 
DATA  SUPPLIED  TO  TASC 


T3486 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

1 

IDENT 

Master  record  type  indicator  (1) 

2-4 

DAY 

Julian  day  (Greenwich  Mean  Time) 

5-6 

HR 

Hours 

7-8 

MIN 

Minutes 

9-10 

SEC 

Seconds 

12 

SITE 

Site  number 

13-15 

DATA ID 

Experiment  number 

17 

SUFSIG 

Master  signal  status:  When  converted  to  a 
binary  number,  first  digit  indicates  whether 
(1)  or  not  (0)  signal  is  present  and  second 
digit  whether  (1)  or  not  (0)  station  is 
blinking 

19 

ENVSKY 

When  converted  to  a  binary  number,  first 
digit  indicates  whether  (0)  or  not  (1) 
envelope  servo  is  enabled  and  second  digit 
whether  (1)  or  not  (0)  skywave  detection  is 
enabled 

20 

SYN INTER 

No  synchronous  interference  (0),  even 
synchronous  harmonic  present  (1),  or  odd 
synchronous  harmonic  present  (2) 

21 

TEMPCOMP 

Temperature  compensation  performed  (1)  or 
not  performed  (0)  in  last  100-sec  interval 

23-27 

DELAY 

Relative  receiver  time  delay  (nsec) 

29-37 

TOA 

Master  time  of  arrival  (0.1  nsec) 

39-42 

RMSTOA 

Master  RMS  time  of  arrival  (nsec) 

44-46 

ETE 

Master  envelope  time  estimate  (0.1  usee) 
(Similar  to  ECD) 

48-53 

SIGLEV 

Master  signal  level  (0.24414  mv) 

56-58 

COUPATTE 

Antenna  coupler  attenuation  (dB) 

60-61 

RECATTEN 

Receiver  attenuation  (dB)  for  master  signal 

63-66 

SNR 

Master  s ignal - to-noise  ratio  (0.1  dB) 
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TABLE  A-4 

AUSTRON-5000  FORMAT  FOR  DATA  SUPPLIED  TO  TASC* 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

2 

STATION 

Station  designator  (M,  X,  Y,  or  Z) 

4-5 

MODE 

Receiver  operating  mode  (AC,  AV,  AF, 
AW,  AS,  AT,  or  K) 

8-10 

GAIN 

Receiver  gain  number 

12-15 

SNR 

S  gnal-to-noise  ratio 

17-21 

CYCLE 

Receiver  cycle  number 

23-30 

TOA 

Time-of-arrival  relative  to  internal 
oscillator  (psec) 

32-39 

TD 

Time  difference  (psec) 

•'Time  is  contained  in  a  separate  record  with  the  following 
format:  DAY(3-5),  HR(7-8),  MIN(lO-ll),  and  SEC(13-14) 
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APPENDIX  B 

FORMATS  FOR  LORAN-C  DATA  DELIVERED  TO  U.S.  COAST  GUARD 


Formats  for  Loran-C  data  delivered  to  the  U.S.  Coast 
Guard  on  half-inch  magnetic  tape  are  presented  in  Tables  B-l 
to  B-5.  The  record  lengths  are  130  characters  for  all  receiver 
formats  and  40  characters  for  the  LPA  format.  All  characters 
which  are  not  included  in  the  tables  are  "blanks." 
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TABLE  B-l 


INTERNAV  LC-204  FORMAT  FOR  DATA  DELIVERED 
TO  U.S.  COAST  GUARD 


—  i  i|i 

1-5 

— 

OBS 

Record  Nuaber 

7 

SITE 

* 

9-13 

RECTYPE 

Receiver  type  (LC204) 

15 

RECNUM 

Receiver  nuaber 

17-19 

DAY 

* 

21-22 

HR 

* 

24-25 

HIM 

★ 

27-28 

SEC 

* 

30-39 

TOl 

Ttae  difference  1  (psec) 

41-46 

DIFCORR1 

Differential  correction  for  TD1  (psec) 

48-57 

TD2 

Tiae  difference  2  (peee) 

59-64 

DIFCORR2 

Differential  correction  for  TD2  (peec) 

66 

IDSHRECD 

* 

68-71 

SMRECDM 

* 

73-76 

SNXECD1 

★ 

78-81 

SNRECD2 

★ 

83 

SWITCH 

* 

85-87 

DATAID 

*■ 

89 

ECSH 

* 

91 

ECS1 

* 

93 

ECS  2 

* 

95 

SUFFSIC 

* 

97-102 

TD1MASJCS 

Editing  Marks  for  TD1 

104-109 

TD2HARXS 

Editing  Marks  for  TD2 

111-118 

TIME 

Tiae  in  sec  since  beginning  of  1979 

* These  data  it cas  ere  described  in  Table  A-l 


TABLE  B-2 

MAGNAVOX  AN/BRN-5  FORMAT  (SECONDARY)  FOR  DATA  DELIVERED 

TO  U.S.  COAST  GUARD 


T-3491 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

1-5 

OSS 

Record  nuober 

7 

SITE 

* 

9-13 

RECTYPE 

Receiver  type  ( BRN5 ) 

15 

RECNUM 

Receiver  nuober 

17-19 

DAY 

★ 

21-22 

HR 

★ 

24-25 

MIN 

* 

27-28 

SEC 

* 

30-39 

TD 

Tiae  difference  1  or  2  (usee) 

41-46 

DIFCORR1 

Differential  correction  for  TD  (nsec) 

48-57 

TOA 

Tiae  of  arrival  1  or  2  (usee) 

59-64 

DIFCORR2 

Differential  correction  for  TOA  (M»ec) 

66-69 

RMSTOA 

* 

71-74 

SNR 

★ 

76-78 

ETE 

♦ 

80-82 

DATAID 

* 

84 

SUFSIG 

★ 

86-91 

SICLEV 

* 

93-94 

RECATTEN 

★ 

96-99 

FRQNT12 

★ 

101 

IDENT 

★ 

103-108 

tdmarjcs 

Editing  narks  for  TD 

110-115 

TOAMARXS 

Editing  aarks  for  TOA 

117-124 

TIME 

Time  in  sec  since  beginning  of  1979 

"These  daca  items  arc  described  in  fable  A- 2 


TABLE  B-3 

MAGNAVOX  AN/BRN-5  FORMAT  (MASTER)  FOR  DATA  DELIVERED 

TO  U.S.  COAST  GUARD 


,49* 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

1-5 

OBS 

Record  Nuaber 

7 

SITE 

* 

9-13 

RECTYPE 

Receiver  Type  (BRN5) 

15 

RECNUM 

Receiver  number 

17-19 

DAY 

* 

21-22 

HR 

* 

24-25 

MIN 

* 

27-28 

SEC 

* 

30-39 

TOA 

Master  time  of  arrival  (wsec) 

4L-44 

RMSTOA 

★ 

46-49 

SNR 

★ 

51-53 

ETE 

* 

55-57 

DATAID 

* 

59 

SUFSIG 

* 

61 

ENVSKY 

* 

63 

SYNINTER 

* 

65 

TEMPCOMP 

* 

67-71 

DELAY 

* 

73-78 

SIGLEV 

* 

80-82 

coupatte 

* 

34-85 

RECATTEN 

* 

87 

ID  ENT 

* 

89-94 

TOAMARKS 

Editing  narks  for  TOA 

96-103 

TIME 

Tine  in  sec  since  beginning  of  1979 

TABLE  B-4 


AUSTRON- 5000  FORMAT  FOR  DATA  DELIVERED 
TO  U.S.  COAST  GUARD 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

1-6 

OBS 

Record  number 

8 

SITE 

Site  number 

10-13 

RECTYPE 

Receiver  type  (AUST) 

15 

RECNUM 

Receiver  number 

17-19 

DAY 

21-22 

HR 

.t 

24-25 

MIN 

* 

27-28 

SEC 

* 

30-37 

TD 

*v 

39-44 

DIFCORR 

Differential  correction  for  TD  (psec) 

46-47 

MODE 

49-51 

GAIN 

* 

53-56 

SNR 

JU 

58-62 

CYCLE 

J. 

64-71 

TOA 

* 

73-78 

TDMARKS 

Editing  marks  for  TD 

80-87 

TIME 

Time  in  sec  since  beginning  of  1979 

TABLE  B-5 

FORMAT  FOR  LPA  DATA  DELIVERED 
TO  U.S.  COAST  GUARD 


CHARACTERS 

DATA  ITEM 

DESCRIPTION 

1-5 

OBS 

Record  number 

7-9 

DAY 

Julian  day  (Greenwich  Mean  Time) 

11-12 

HR 

Hours 

14-15 

MIN 

Minutes 

17-18 

SEC 

Seconds 

20-23 

LPA 

Local  phase  adjustment  (nsec) 

25 

STATION 

Loran-C  transmitter  (X,  Y,  Z) , 
or  "B"  or  "E"  to  indicate 
beginning  or  end  of  segment  of 

LPAs  for  which  time  is  uncertain 

27 

MASTER 

Set  to  if  LPA  is  applied  through 

master  and  other  two  secondaries 

29-36 

TIME 

Time  in  sec  since  beginning  of  1979 

*The  LPA  times  are  known  to  lie  between  the  B-time  and  E-time. 
LPA  times  for  the  records  between  the  B-  and  E-records  are 
set  to  the  B-time. 
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APPENDIX  C 

PARTIAL  LISTING  OF  LORAN-C  DATA  BASE 
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Examples  of  the  data  formats  described  in  Appendix  B  ^ 

are  presented  in  Figs.  C-l  to  C-5.  Lead  and  final  file  records 
are  not  included  in  these  listings,  but  are  included  on  the  j 

magnetic  tapes.  Column  labels  are  provided  here  for  convenience,  j 

but  are  not  included  on  the  magnetic  tapes. 
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Figure  C-l  Partial  Listing  of  Internav  LC-204  Data 
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Figure  C-4  Partial  Listing  of  Austron-5000  Data 
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